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We have modelled the high-energy astrophysics of the inner 200 pc of the Galaxy with a view to 
,_i explaining the diffuse, broad-band (radio continuum to TeV 7-ray), non-thermal signal detected 

from this region. Our modelling pins down the ISM parameters for the environment wherein 
X cosmic ray (CR) electrons and ions reside in the Galactic centre (GC). We find that the magnetic 

(-H field in this region is 100-300 /zG, the gas density < 60 cm^^, and that a powerful (> 200 km/s) 

0^ 'super'-wind acts to remove > 95% of the cosmic rays accelerated in the region before they have 

Q time to lose their energy in situ. The ^ 10^^ erg/s carried away by the GC cosmic ray protons is 

precisely enough to energise the ~GeV y-ray emission from the Fermi 'bubbles' recently found to 
^ extend north and south of the GC out to distances of ~^10 kpc, provided that the bubbles constitute 

thick targets to the GC protons and that the situation has reached steady state. In such a situation 
of 'saturation' the hard, uniform spectrum of the bubbles are explained and secondary electron 
^ synchrotron explains the non-thermal microwave emission found in WMAP data mirroring the 

bubbles. Given the very low density of the bubble plasma (<0.01 cm^^*), the pp loss time in the 
Bubbles is > 5 Gyr Our scenario thus has the startling implication that a GC source of non- 
thermal particles of time-averaged power 10^^ erg/s has persisted since the youth of the Galaxy. 
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1. Introduction 

The Galactic centre (GC), besides the intrinsic interest it holds, provides an interesting, poten- 
tial analogue to the nucleus of a luminous star-burst galaxy. Indeed, the ISM conditions prevailing 
in the inner ~ 100 pc (in radius) of the Galaxy - the region under consideration here - render it 
arguably more akin to the environs of a star-burst than than to the relatively quiescent conditions 
of the Galactic disk. In particular, the GC contains something hke 5% of the Galaxy's molecular 
hydrogen allocation [19] implying a very high, volumetric-average gas density in the region. More- 
over, the energy-densities of the various GC ISM components are 1-2 orders of magnitude larger 
than those found locally, as is the areal density of star-formation and attendant supernova activity. 
For instance, as we have recently shown [5], the GC is threaded by a remarkably strong magnetic 
field of ~100 /iG (cf. with ~5 liG for the Galactic disk). Such a field imphes that >10% of the 
Galaxy's magnetic energy is contained in only <0.05% of its volume. 

In recent years a picture has begun to emerge that the inner regions of star-forming galaxies 
should i) be important sources of 7-rays in the universe [21, 28, 12]; ii) drive powerful galactic 
winds [12] and iii) therefore, be important shapers of the inter-galactic medium, particularly its 
metallicity [25]. Interestingly, recent observations reveal that GC is also a significant 7-ray source 
exhibiting both point-like GeV [3] and TeV [1] emission coincident with Sagittarius A* (at the 
dynamical centre of the Galaxy and most likely associated with the super-massive black hole found 
there) and diffuse emission also at both GeV [9] and TeV [2] energies. 

We argue here for a further similarity to star-bursts: many direct observations (see Appendix 
C of [6] and references therein) and, in addition, the non-thermal evidence reviewed below point 
to the existence of a power outflow or wind of at least a few 100 km/s out of the GC. We finally 
consider the implications of the injection of this wind fluid - composed of very hot plasma, cosmic 
rays and, presumably, 'frozen-in' magnetic field lines - into the Galaxy-at-large and explore a 
compelling connection between this outflow and the recently-discovered [26, 1 1] 'Fermi Bubbles'. 

2. Broad-band modelling of the Galactic centre 

We have created a one-zone model of the injection, cooling, and escape of relativistic pro- 
tons, electrons, and secondary electrons (and positrons) from the inner ~100 pc (in radius) of the 
Galaxy. This is the approximate region for which the HESS telescope has reported [2] a diffuse, 
TeV 7-ray flux. To match the TeV data we have collected archival ~GHz radio continuum data 
covering the same region (see Appendix D of [7] for radio data sources) and also GeV data [18]. 
Unfortunately, the latter is heavily polluted by the contribution of i) individual point sources within 
the field [4] including a GeV source positionally coincident with Sgr A and the GC TeV point 
source [3] and ii) by emission from CRs in the line-of-sight along the Galactic plane but out of 
the GC [29] and does not usefully constrain our broadband modelling of the diffuse emission from 
the region. As justified elsewhere [7], in our modelling we assume that the particle astrophysics 
can be accurately described to be in quasi-steady state and that the particle transport timescale is 
energy-independent. We also assume homogeneity and isotropy in our modelling. 

Our modelling approach is to find - as a function of environmental and other parameters - the 
steady-state populations of relativistic protons and electrons within the region of interest. (Note 
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that in our modelling we neglect for simplicity the poorly-constrained ionic component of the 
CR hadronic population heavier than protons.) We then self-consistently determine (for the same 
environmental parameters) the radiative output of these populations. Relevant radiative processes 
are neutral meson decay for protons and synchrotron, inverse-Compton and bremsstrahlung for 
electrons. We self consistently track both primary and secondary electron emission in our radiative 
modelling. Finally, we use a minimization procedure to determine the parameters describing 
the proton and electron populations whose radiative outputs in the given ISM environment give the 
best fit to the ~ GHz radio continuum spectrum and the ~TeV /-ray spectrum detected from the 
HESS field (as previously noted, the GeV data only define upper limits to the broadband emission 
from the region and we do not attempt to fit them directly). 

Environmental parameters that vary within our modelling are magnetic field B, ambient hy- 
drogen number density nu (in whatever form), and the energy-independent timescale over which 
particles are advected from the system, tesc = h/v^nd, where /j = 8 kpc tan(0.3°) ~ 40 pc. (For the 
reasons explained at length in ref. [7] advection must be the dominant CR transport mechanism in 
the GC environment so we do not consider CR diffusion in our modelling) 

We assume that protons and electrons are injected into the GC ISM with distributions governed 
by power laws in momentum (with identical spectral indices jp = = 7)- The relative normaliza- 
tion of the injection distribution of electrons to that of protons (at relativistic energies) is given by 
the coefficient K^p which is also left as a free parameter in our modelling. The absolute normaliza- 
tion of the distribution of protons at injection is, finally, specified by Qp (in units eV^^cm^^ s^'), 
also a free parameter within the model. 

3. Implications of Galactic centre modelling 

Our fitting procedure finds acceptable fits for magnetic field amplitudes around ~ 100 
wind speeds around v„ind ~ few x 100 km/s, total power in all non-thermal particles of ~ 10^^ erg/s, 
gas densities ai^ound nu ~ 10 cm^^, injection spectral indices y ^lA, ionization rates of ~ 10^'^ 
s^^ and electron to proton ratios (at injection at 1 TeV) of K^p ~ 10^^. Below we illustrate (figs. 1 
the modelled broadband spectrum for the best-fit case (which achieves — 7.9 for dof = 9). 

On the basis of our modelling we can determine a number of interesting facts about the GC. 
Firstly, the totality of non-thermal signals require a contribution from both primary electrons and 
protons. Neither scenarios where primary electrons alone (in which the observed TeV emission is 
mostly provided by IC emission) nor where primary protons alone (in which secondary electrons 
supply the observed synchrotron radiation) provide acceptable fits to the data. This latter may be in 
contrast to the case presented by luminous star-bursts where it has been claimed [27] that secondary 
electrons probably do supply most of the observed synchrotron emission. 

Secondly the gas environment where the non-thermal radiation is being generated is less dense 
than ~ 60 cm^^ (at 2a confidence) with a best fit value close to 1 cm^^. Given that the upper end 
of the allowed nu range is less than the volumetric average gas density through the region (~ 120 
cm^^) this is an indication that CRs - even the > 10 TeV protons responsible for generating the 
TeV 7-ray emission do not penetrate into the densest gas in the region (where star-formation is 
occurring). This is, again, apparently in contrast to the situation presented by star-bursts where, it 
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Figure 1: HESS field broadband spectrum energy distribution for the best-fit case. Fitted parameters include 
magnetic field amplitude 200 /iG), ambient gas density (6 cm^^), at-injection electron-to-proton ratio at 
1 TeV (~0.005), and wind velocity (~400 km/s). Curves are divided into (i) dashed - primary electron 
emission; (ii) dotted - secondary electron (and positron) emission; and (iii) solid - total emission at a 
given photon energy. Emission processes are: blue - synchrotron; red - bremsstrahlung; green - inverse 
Compton; and brown, dot-dashed - neutral meson decay. The total 7-ray flux is shown in purple. Data 
are from (at Ey - lO^^) radio, (at £y - 2 keV) Ginga [30] with , (at Ey ~ 1 GeV = l(f eV) Fermi [18], 
and (at Ey^ I TeV = lO'^ eV) HESS observations [2]. We only display modelled synchrotron emission at 
radio wavelengths; other processes that combine with this to give the region's total radio emission are not 
pictured. Finally, note that the Fermi data points define only upper limits to the diffuse emission. 



has recently been claimed [20], CRs modify the gas conditions and chemistry where star-formation 
is occurring, potentially biasing the initial stellar mass function towards more massive stars. 

This result - that GC cosmic rays do not penetrate into the dense gas - is consistent with 
another finding: a powerful, star-formation driven superwind blows out of the region with a speed 
of > 200 km/s (see §3.1 below). Given the wind, comparison of relevant timescales indicates that 
CRs do not remain long enough in the region to penetrate into the dense molecular gas. 

We note that the CRs, however, do apparently constitute important sources of heat and ion- 
ization for the warm, diffuse molecular gas phase enveloping the molecular gas cores in the GC 
environment. 

3.1 Existence of an outflow 

Even without detailed modelling, evidence for a strong outflow can be gathered from a com- 
parison of the thermal and non-thermal signals detected from the inner ~ 100 pc (in radius) of the 
Galaxy [8]. At TeV energies, the luminosity expected from the inner 1.6° of the Galaxy were the 
system to be calorimetric to all (sufficiently energetic) cosmic rays is given by: 

- i X X /„ X 2 X 10" erg/s .co ( j,;,^,^) (^) O- D 

where tJcr is the fraction of supernova mechanical power that goes into non-thermal protons nor- 
malized to T]o.io = T]cr/0.1 and /xev — 0.05 is the fraction of total CR proton power in protons 
sufficiently energetic to generate TeV 7-rays (adopting ~ 10 TeV for the mean energy of the parent 
proton of a TeV 7-ray and spectral index 7= 2.3 in agreement with the TeV spectrum). The ground- 
based 7-ray telescope HESS has detected [2] hard spectrum, diffuse ~TeV 7-rays from precisely 
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the same inner 1.5° region as noted above. This flux, though interpreted originally in the context 
of an explosive injection of cosmic ray hadrons [2] in a single GC event a few thousand years ago, 
nevertheless proffers an upper limit on the steady-state /-ray intensity from the region. The total, 
>TeV luminosity of the region is 1.2 x 10^^ erg/s, which, given the above, is only about 1% of the 
luminosity expected were the system calorimetric for cosmic ray protons. Thus the vast majority 
of the power injected into GC cosmic rays must be carried outside the region. The same conclusion 
can be reached from consideration of the non-thermal radio continuum emission from the region [8] 
and even from the heavily-polluted GeV emission; both these demonstrate that at least 90% of the 
power injected into the region's non-thermal proton and electron populations is lost non-radiatively 
from the region. The GC, then, loses r\cR x "^sn x Esn ~ 0.1 x 0.04/century x 10^°erg ~ 10^^ erg/s 
in hard-spectrum CRs into the general Galactic environment. 

At the heuristic level presented here the conclusion that there is a strong wind out of the region 
would seem to require that i) we have correctly determined the supernova rate in the region and 
that, ii) indeed, lO^*' erg per supernova is injected into non-thermal particle populations. What 
if our supernova rate is off-beam or, for whatever reason, GC supernova remnants (SNRs) are 
comparatively inefficient CR accelerators? Would this not, then, imply that we could no longer 
infer the existence of a wind? In fact our detailed modelling generates an independent constraint 
on the power injected into CRs in the region. At > 10^^ erg/s, fortuitously or not, this power is 
essentially what we expect given our supernova rate determination with lO^*' erg per supernova 
injected into non-thermal particles assumed. So the conclusion that there is a star-formation-driven 
wind out of the region seems to be robust. What is more, we reach also reach the interesting 
conclusion that GC SNRs are at least as efficient as the Disk variety as CR accelerations. 

4. Connection to the 'Fermi bubbles' 

Recently NASA announced the startling discovery [11, 26] by Fermi of two enormous gamma- 
ray emission structures that hang like lightglobes above and below the centre of the Milky Way. 
These 'Fermi bubbles' extend an astounding 10 kpc from the plane of the Galaxy. At lower Galactic 
latitudes these structures are coincident with a non-thermal microwave 'haze' found in WMAP 20- 
60 GHz data [14, 10] and an extended region of diffuse X-ray emission detected by ROSAT [23]. 
Thus far the Bubbles have been typically understood as illuminated by a mysterious population of 
youthful and highly energetic (~TeV) electrons of age 10 million years or so which simultaneously 
inverse-Compton radiate (off the CMB) at ~GeV energies and synchrotron radiate at microwave 
frequencies. However, given the severe radiative energy losses experienced by electrons, the hard 
spectrum, uniform intensity, vast extension, and energetics of the bubbles render the origin of this 
particle population extremely mysterious [14, 10, 11, 17, 26]. In particular, even accounting only 
for energy losses on the CMB, transport of >TeV, IC -radiating electrons to the requisite distances 
from the plane would require velocities of > 0.03 c, too fast for a Galactic wind (though an AGN 
jet potentially offers a suitable delivery mechanism in this connection [15]). 

A viable alternative to the idea that CR electrons generate the Bubble /-rays is, however, that 
they originate in the collisions between cosmic ray ions and the Bubbles' low-density plasma. 
In fact, the wind out of the GC we have identified above canies a power of ~ 10^^ erg/s out of 
the region, as we have emphasised, and this is precisely enough, in steady state and assuming 
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the Bubbles represent thick targets to the injected protons, to sustain the observed 1-100 GeV 
luminosity of the Bubbles of 4 x 10^' erg/s. This scenario explains many aspects of the Bubbles' 
non-thermal phenomenology. Firstly, the hard spectrum of the /-rays is explained: in contrast to the 
situation in the Galactic disk where energy-dependent confinement implies a steepening of the in 
situ spectrum of CRs away from their injection distribution, the Bubble CRs, trapped independently 
of energy by hypothesis, follow their injection distribution. Likewise the hard-spectrum 'WMAP 
haze' [14, 10], coincident at lower Galactic latitudes with the ~ GeV y-ray emission, is explained 
in this scenario as a result of secondary electron synchrotron emission which, again, would provide 
a signal of precisely the right luminosity. (We note in passing that decay of charged mesons - 
which leads to secondary electrons and positrons - also produces high-energy neutrinos and that 
- if our scenario is correct - the Bubbles should also be a significant source for a future, km^ 
class. Northern Hemisphere neutrino telescope.) Finally, a robustly-detected [26] down-turn in the 
Bubbles' SED below ~GeV is explained as a natural result of the kinematics of neutral pion decay. 

The 'cost' of this scenario is the very long timescales implied: given the low-density of the 
Bubble plasma, < 0.01 cm^"*, the pp loss time blows out to > 5 Gyr and the Bubbles are required 
to have existed for at least this time (in order that steady state be satisfied) and to effectively trap 
CR protons (up to at least ~TeV) over the same timescale. Thus one requires that the GC has 
sustained injection of ~ 10^^ erg/s in CRs into the base of the Bubbles for multi Gyr timescales. 
Such might seem hard to credit but is not unreasonable upon further reflection: the morphology of 
the Bubbles cleai^ly privileges the GC and the GC is perhaps the single, spatially-localized site in 
the Galaxy where SF over multi-Gyr timescales is assured [22, 13]. Note that while by no means 
required by our scenario, it is interesting that the current level of GC star-formation - and resulting 
cosmic ray luminosity - is close to the time-averaged value required in our scenario for the origin 
of the Bubbles. This suggests a system in steady state. 

A final couple of notes are in order. Firstly, aside from the sustained star-formation occurring 
in the GC, an a priori suspect to ultimately power the Bubble emission is the central, supermassive 
black hole (SMBH) [26, 15]. In the sort of hadronic scenario we have explored this would - 
just as for star-formation and concomitant supernova activity - be required to generate a time- 
averaged power of ~ 10^^ erg/s in CRs. It is interesting in this context, then, that such a cosmic ray 
luminosity is rather close to the minimum required by analysis of the central, point-like GeV 7-ray 
source under the assumption that it is hadronic in nature [3]. 

Secondly, we note that Socrates et al. have recently postulated [24] the existence of an 'Ed- 
dington limit in cosmic rays' following the logic that, as cosmic rays diffuse outwards from some 
central source, they exchange momentum with surrounding gas via scattering on 'frozen-in' mag- 
netic field inhomogeneities and may, therefore, arrest gas accretion beyond some limiting luminos- 
ity. This is roughly estimated to be Z.^^^ ~ lO^^LEdd- For the case of the GC SMBH (of mass 
~ 4 X 10^ M0) or the central nuclear star cluster (of mass ~ 3 x 10^ M©: [16]) the GC's Central 
Massive Object has a rough 'cosmic ray Eddington limit' in the range (0.5 — 4) x 10^^ erg/s, again 
interestingly close to the value required to sustain the Bubbles' current 7-ray emission. In conclu- 
sion, over the last few billions years, the GC seems to have emitted a time-averaged cosmic ray 
luminosity of close to the maximally- allowed value consistent with sustained accretion. We finally 
remark that, regai^dless of whether it is GC stai^-formation or low-level but sustained activity of the 
SMBH that energises the Bubbles, in our hadronic scenario these remai^kable structures constitute 
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a perfect calorimetric recording of Galactic centi-e activity over the history of the Milky Way. 
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